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TECHNICAL NOTE 322? 


APPLICATION OF TWO-DIMENSIONAL VORTEX THEORY TO THE 
PREDICTION OF FLOW FIELDS BEHIND WINGS OF WING -BODY 
COMBINATIONS AT SUBSONIC AND SUPERSONIC SPEEDS 
By Arthur Wm. Rogers 

SUMMARY 


A theoretical Investigation has heen made of a general method for 
predicting the flow field Behind the wings of plane and cruciform wing 
and body combinations at transonic or supersonic speeds and slender con- 
figurations at subsonic speeds. The wing trailing-vortex wake is repre- 
sented initially by line vortices distributed to approximate the spanwise 
distribution of circulation along the trailing edge of the exposed wing 
panels. The afterbody is represented by corresponding image vortices 
within the body. Two-dimensional line-vortex theory is then used to 
compute the induced velocities at each vortex and the resulting displace- 
ment of each vortex is determined by means of a numerical stepwise inte- 
gration procedure. The method was applied to the calculation of the 
position of the vortex wake and the estimation of downwash at chosen tail 
locations behind triangular -wing and cylindrical-body combinations at 
supersonic speeds. The effects of such geometric parameters as aspect 
ratio, angle of attack and incidence, ratio of body radius to wing semi- 
span, and angle of bank on the vortex wake behind wings of wing -body 
combinations were studied. The relative importance of wing vortices, 
the corresponding image vortices within the body, and body crossflow in 
determining the total downwash was assessed at a possible tail location. 

It was found that the line -vortex method of this report permitted 
the calculation of vortex paths behind wings of wing-body combinations 
with reasonable facility and accuracy. A calculated sample wake shape 
agreed qualitatively with one observed experimentally, and sample results 
of the line -vortex method compared well with an available exact crossflow- 
plane solution. An empirical formula was derived to estimate the niomber 
of vortices required per wing panel for a satisfactory computation of 
downwash at tail locations . It was found that the shape of the vortex 
wake and the ultimate number of rolled-up vortices behind a wing depend 
on the circulation distribution along the wing trailing edge. For the 
low-aspect-ratio plane wing and body combinations considered, it appeared 
that do^mwash at horizontal tail locations is largely determined except 
near the tail-body juncture by the wing vortices alone for small ratios 
of body radius to wing semispan, and by the body upwash alone for large 
values of that ratio. 
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INTRODUCTION 


Satisfactory aerodynamic design of high-speed aircraft requires 
knowledge of the interference flow field resulting from wing -body- tail 
interaction. The behavior of the wing vortex wake in the presence of 
the body directly affects the air stream flowing past the tail svirfaces. 
In particular, for certain relative sizes and positions of wing, body, 
and tail, the wing-body vortex wake produces large stream angles at the 
tail surfaces and nonlinear variation with angle of attack of the pitch- 
ing moment contributed by the tail. Therefore, it is desirable to be 
able to predict flow fields at the tail location for a given configura- 
tion in order to evaluate stability and control requirements. 

Furthermore, vapor-screen studies in supersonic wind tunnels have 
shown that the configuration of the vortex wake from the wings of some 
wing-body combinations differs markedly from the usual conception of a 
sheet whose side edges curl over to form a single pair of vortices. 

For example, fi^re 1(a) shows the S-shaped cross section of a vortex 
sheet from each panel of a high-aspect-ratio supersonic wing and body 
combination at angle of attack.^ Here the ultimate vortex pattern far 
downstream of the wing consists of two vortices from each wing panel, 
both rotating in the same direction. Such a vortex wake results in a 
flow field different from the more well-known pattern, and the load on 
a tail situated in that flow field differs from the tail load associated 
with a single-vortex wake. It is important, therefore, to determine 
the conditions under which this unusual wake pattern occurs. 

Considerable investigation (refs. 2-12) has been devoted to the pre- 
diction of flow fields behind wings, both plane and cruciform. Refer- 
ence 2 uses linearized conical-flow theory to calculate sidewash and 
downwash for a flat vortex sheet at the wing trailing edge and in the 
Trefftz plane, references 3 and k use supersonic potential doublet dis- 
tributions to determine downwash in the plane of a flat wake and in the 
vertical plane of symmetry for any distance downstream of a wing. Ref- 
erence 5 employs pressure doublets to obtain general expressions for 
the induced velocities in space behind a wing, and references 6 through 8 
utilize line -vortex theory to predict sidewash and downwash in the flat- 
vortex wake and in space. The flow in the Trefftz plane behind super- 
sonic wings is treated in reference 9 ly means of linearized conical-flow 
theory. The downwash based on a flat deflected vortex sheet and the 
downwash based on two deflected vortex lines are calculated in refer- 
ence 10. The rolling up of the trailing- vortex sheet behind wings is 
analyzed in reference 11, and reference 12 considers the motion of single 
vortices from each panel of a cruciform wing. The literature on the 
vortex wake behind wings of wing -body combinations, however, is still 
relatively meager (refs. 13 and l4) . For wing-body combinations, ref- 
erence 13 accounts for the effect of the fuselage on the flow field by 


^A description of the vapor-screen technique is given in reference 1. 
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considering the motion of two fully rolled-up vortices in the presence 
of a circiilar cylinder. Reference l 4 treats a specific configuration by 
slender-body theory with the two assumptions of either a flat vortex 
sheet or two fully rolled-up vortices. 

The general problem to be considered here is the determination of 
the interference flow field behind wings of wing-body combinations 

fig* 2) at subsonic or supersonic speeds. The dual purposes of 
this report are to present in detail and evaluate the method outlined 
in reference I3, and to apply it to the calculation of vortex paths and 
to the effects of wing aspect ratio, span loading, angle of attack, ratio 
of body diameter to wing span, angle of bank, and wing incidence on vor- 
tex paths and downwash at a possible tail location. 

Part I of this report presents the results of the application of 
the method. All calculations of span loading were made for triangular- 
wing and cylindrical-body combinations at a Mach nvunber of 2 . 0 . Wing 
aspect ratios of 2/3, 2, and if are considered, combined with bodies 
whose radii are 0 . 2 , O.h, 0 . 6 , and 0.72 times the wing semispan. The 

selected for downwash computations in this part represents 
a missile-type configuration, for which the horizontal tail is assumed 
to be located in the body diametral plane 10 body radii downstream of 
the wing trailing edge. 


Part II of this report contains a detailed description of the line- 
vortex method. An examination is made of the effect on downwash at 
chosen tail locations of the ntmiber of vortices used to represent the 
trailing vorticlty. The nature of the stepwise integration method used 
is discussed, and solutions obtained by the stepwise integration method 
are compared with known exact solutions. Downwash is also computed at 
the tail of a high-tall airplane-type configuration, characterized by a 
horizontal tail 3 body radii downstream of the wing trailing edge and 
2 body radii above the body axis. 


SYMBOLS 


A aspect ratio 

c local wing chord 

wing root chord in plane of syrmnetry 

directed line segment of a contour 
E elliptic integral of second kind 

unit vectors in x,y,z directions 


k 
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k 

L 

M 

n 

P 

Ap 

q 

q 

r 

s 

Sm, 

s* 

t 

Pm 

u,v,w 

Vo 


aspect-ratio correction factor defined in Appendix B 
lift 


Mach mmber 

nvunber of vortices in summation 
static pressure 

pressiare differential across wing surface, p^ - p^ 


free-stream dynamic pressure. 



velocity vector, uT + vj + wk 
Body radius 

local horizontal-wing semispan 

maximum horizontal -wing semispan 

spanwise position of wing vortex 

local vertical-wing semispan 

maximxmi vertical -wing semispan 

velocities in x,y,z directions (See fig. 2.) 

free-stream velocity 

chordwlse distance from leading edge 

Cartesian coordinates of wind axes (See fig. 2.) 


_ load coefficient 

q 

a angle of attack of "body axis, radians 

3 - 1 

3’ angle of yaw, a sin •?, radians 

r(y) spanwise distribution of circulation (defined by eq. (Al)) 

A increment 
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e 

downvash angle, — ^ 
Vo 

e 

tan-^ f 

A 

ving leading -edge sveep angle ^ 90 ^ - w, radians 


dimensionless rectangular coordinates of wind axes^ 
(See fig. 2.) 

P 

air density of free stream 

a,T 

coordinates of source -point in xy plane 

9 

angle of lank^ radians 

• • ' 

► velocity potentials 

(a) 

wing semiapex angle, tan“^ radians 

Cr 

Subscripts 

B 

body 

j 

indices of summation 

1 

lower 

m 

maximum 

u 

upper 

W 

wing 

LT 

linear theory 

SBT 

slender-body theory 
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I.- APPLICATION OF METHOD TO THE PREDICTION OF VORTEX PATHS 
AND DOWNWASH BEHIND TRIANGULAR -WING AND CYLINDRICAL-BODY 
COMBINATIONS AT SUPERSONIC SPEED 


General Description of Method 


The method used herein follows that of reference I3. The circula- 
tion distribution along the wing trailing edge of a given wing-body con- 
figuration is required as the initial condition. In this report, the 
circulation distribution is determined by modified slender-body theory. 
The circulation distribution is then replaced by a finite number of vor- 
tices which trail downstream and represent the wake vorticity. The 

of" the afterbody is accounted for by vortices placed within the 
body at the image position of each of the wake vortices. It is assumed, 
in general, that the flow changes in the stream direction are such that 
in the velocity potential equation 


(1 - M^) + 9yy + - 0 

the first term is negligibly small, that is, 

(1 - M^) cpxx = 0 

so that 


Cp + CD 

yy ^ 


zz 


= 0 


Thus the solution is independent of Mach number. Such an assumption 
of two -dimensionality is valid for slender, pointed wings and bodies at 
subsonic speed, and at supersonic speed when the entire plan form lies 
within the body nose Mach cone, or near M = 1 for more general con- 
figurations.^ 


Consistent with this slender-body theory assumption, the bound vor- 
tices within the wing are neglected in computing vortex paths and down- 
wash. The trailing vortices are extended to infinity upstream and down- 
stream, and the induced velocities in crossflow planes are calculated by 
two-dimensional line -vortex theory and body potential crossflow. In the 
crossflow plane, therefore, the analysis considers the following picture: 


^See reference 15 for a discussion of these limitations. 
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The motion of each of the vortices in the crossflow plane is computed, 
and the results applied to the wing-hody problem by relating time in the 
crossflow plane motion to distance downstream of the wing trailing edge. 

■nie wing and image trailing vortex system used in the stepwise cal- 
culations constitutes the only vortex pattern considered here. It is 
known that at high angles of attack and for large ratios of body diameter 
to wing span, viscous crossflow produces vortices above the body (e.g., 
fig. 1 and ref. l) which significantly affect the flow field. A thorough 
understanding of the mechanism of formation, strength, and stability of 
crossflow vortices is as yet lacking. Consequently, no attempt has been 
made to account for them. 

On the basis of a large number of computations made during this 
investigation, the following analysis was made of the effect of geometric 
parameters on the vortex wake and the induced flow field behind wings of 
wing-body combinations . Although only triangular wing and body combina- 
tions were studied, the general conclusions should be applicable to any 
wing plan form for which the circulation distribution is similar to the 
types contained herein. 
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Effect of Reduced Aspect Ratio, PA, and Span Loading 
on the Rolling Up of the Vortex Wake 


Figure 3 shows the calculated detailed rolling up of 10 vortices 
representing a typical wake behind a low-aspect-ratio configuration. 
Vortex coordinates for this figure are listed in table I. It can be 
seen that at 10 body radii downstream of the wing trailing edge, 7 of 
the 10 vortices have already spiraled into a single group which trails 
downstream in approximately the free-stream direction. The circulation 
distribution along the trailing edge of the subsonic-leading-edge wing 
of this configuration is of almost elliptic shape. The vortex sheet 
behind a wing of a high-aspect-ratio combination, represented in fig- 
ure 4 , distorts much slower toward its final configuration and shows a 
different pattern. Vortex coordinates for this figure are listed in 
table II. The different shape of the vortex sheet is due to the change 
in span loading with reduced aspect ratio PA. 

The difference in span loading for the high- and low-aspect-ratio 
configuration is a consequence of the supersonic or subsonic leading 
edge, respectively. As a consequence of its supersonic leading edge, 
the wing trailing edge of the high-aspect-ratio combination of figure 4 
has a theoretical circulation distribution which is linear from about 
the mid-semispan to the tip. This linear distribution is represented 
by the eight equal, uniformly spaced vortices shown, the remaining two 
vortices arising from the slender-body type of loading (eq. (3) of 
part II ) inboard of the intersection of leading-edge Mach line and the 
wing trailing edge. Now, it is characteristic of a uniform distribution 
of isolated vortices along a line, such as those shown in figure 5^ that 
the rolling up proceeds in the form of a symmetrical, S-shaped curve 
rotating about the centroid. The final configuration of the vortices in 
figure 5 will consist of two equal vortex cores, each containing four 
vortices, rotating symmetrically about the fixed centroid. In figure 4 , 
therefore, the S-shaped rolling up occiars for the uniform portion of the 
sheet, although. the symmetry about and fixity of the centroid in space 
does not occur because the flow field is due not only to the 8 vortices 
but also to the other 12 wing vortices, 20 image vortices, and the body 
crossflow. The inboard portion of the vortex sheet displaces downward 
and outward in the conventional manner of figure 3* In summary, then, 
the effect of aspect ratio and the consequent change in circulation dis- 
tribution is to change the shape of the rolling-up vortex sheet as in 
sketch (b) . The final vortex pattern for the high-aspect-ratio case 
depends on the extent of the span over which the circulation distribu- 
tion along the wing trailing edge is linear. A comparison of the calcu- 
lated S-shaped sheet and an experimentally observed vortex sheet is 
shown in figures 1(a) and 1(b). 
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Low 3A HigL 3A Very high 3A 

Trailing-edge circulation distribution 



Shape of rolling-up vortex sheet 



Final vortex pattern 
Sketch (h) 


Effects of Geometric Parameters on Vortex Paths 
and Downwash at a Tail Location 


In the following sections, calculated vortex paths and downwash are 
presented for a wing -body -tail combination at a Mach n\imber of 2.0. The 
tail location at which downwash is determined is 10 body radii downstream 
of the wing trailing edge, a possible missile -type configuration. The 
vorticity shed from each wing panel is represented by a single vortex. 

For the comparison of the relative importance on downwash at the tail of 
the wing vortex, image vortex, and body crossflow, a single vortex for 
each wing panel suffices, in general, although the computation of the 
actual magnitude of doirawash requires several vortices. In part II of 
this report an investigation is made of the effect of the n\Mber of vor- 
tices used to represent the trailing vorticity on the computed downwash 
at chosen tail locations. 

Effect of aspect ratio .- Figures 6(a), (b) , and (c) present the 
effect of aspect ratio on downwash, calculated at a missile-type tail 
location, for a small body and wing combination (r/sm = 0.2) at an angle 
of attack of 5°. For the low-aspect-ratio combination, the approximation 
that the total downwash is caused by wing vortices alone is good except 
in close proximity to the tail-body juncture. For increasing aspect 
ratios the approximation is good only at greater spanwise distances from 
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the body. The reason for this result is simply that the downwash from 
the wing vortices decreases with increasing aspect ratio for the condi- 
tion of a given ratio of hody radius to wing semispan. The wing lift 
decreases with increasing aspect ratio because of span load changes, and 
the effect of the hody upwash therefore becomes more pronounced. Slender- 
body theory predicts that the wake downwash angle exactly cancels the 
flow angle of attack inboard of the tip. It is interesting to note that 
inboard of the vortex location for the aspect ratio 2/3 wing (A ss 4) the 
average downwash angle is roughly the negative of the angle of attack of 
the wing -body configuration. 

Effect of angle of attack .- Vortex paths behind wings of two 
triangular-wing and body combinations at various angles of attack are 
presented in figure 7. It is seen that increased angle of attack results 
in a more pronounced inward and downward motion of the wing vortex with 
increasing distance p downstream, although for low-aspect-ratio com- 
binations the initial vertical motion is upward due to body upwash. 

Farther downstream the distance between the body and the vortex wake 
increases, with a corresponding decrease in the effect of body upwash. 

The vortices from the wing panels of the low-aspect-ratio configuration 
then move downward more rapidly because of their nearness to each other. 

Examination of the equations of motion of the vortex wake, presented 
in part II of this report, shows that for a given wing -body combination 
a single set of computations can be made for all angles of attack a. 

The results of figures 7(a) and (b) are replotted against pa in fig- 
ure 7(c), which then applies to all angles of attack, at the same Mach 
number . 

The effect of angle of attack on downwash is illustrated by fig- 
ures 8(a), (b), and (c) for an aspect ratio 2/3 wing and large-body com- 
bination, and by figures 8(d), (e), and (f) for an aspect ratio h wing 
and small -body combination. For the low-aspect-ratio configuration, it 
appears that the largest portion of the total downwash is contributed by 
the body upwash because of the large body. The difference in the shape 
of the total downwash c\irves of figures 8(a), (b), and (c) is due to the 
fact that the vortex wake passes farther above the horizontal tail as 
the angle of attack increases. For the high-aspect-ratio configuration 
the total downwash is not primarily caused by any single component 
because of the small body size. Wing and image vortices together with 
body upwash must be considered in calculating total downwash at any 
angle of attack for such configurations. 

Effect of ratio of body radius to wing semispan .- Another parameter 
governing the behavior of the vortex wake is the ratio of the body radius 
to the wing semispan, r/sm- Figure 9 presents the results for an aspect 
ratio 2/3 wing of fixed span, alone and in combination with bodies of 
‘iifi’si’snt radii. The variation of the vortex strength with body size, 
given by equation (3) with y = r, is plotted in figure 10. The vortex 
strength goes to zero for r/sm = 1 since the entire wing then is 
enclosed by the body. It is seen in figure 9 that the initial slope of 
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the vertical displacement of the vortex path near the wing trailing edge 
increases in magnitude as r/sm increases. This is the result of 
increasing upwash in the body crossflow field coupled with correspond- 
ingly decreasing wing vortex strength (see fig. lO) . The final uniform 
downward motion of the vortex pair begins when the body is far removed 
from the vortices, and the path lies above that for the wing alone case 
in which the vortex pair moves linearly downward from the trailing edge, 
^e lateral motion of the trailing vortex is shown in the lower part of 
igure 9* Together with the decreased strength of the wing vortex for 
large values of r/sm, there is an initial outboard shift of the vortex 
at the trailing edge. For the isolated elliptically loaded wing, y/sjri 
equals rt/4 and it approaches 1 (the initial vortex location moves 
toward the wing tip) as r/sm increases from 0 to 1. The lateral vor- 
tex motion downstream is affected by the sidewash component of the body 
crossflow field. The larger the body for a given wing, the more rapid 
IS the inboard motion of the vortex. 

Figures 6(a), 11(a), 8(a), and 11(b) show the effect of r/sm on 
downwash for ratios of body radius to wing semispan of 0.2, 0.4, 0.6 
and 0.72, respectively. The results shown in the figirres can be antici- 
pated qualitatively. For an r/sm of 0.2, figure 6(a) shows that the 
downwash contributed by the image vortices almost entirely cancels the 
body upwash, except in close proximity to the body- tail juncture. Thus, 
for small ratios of body radius to wing semispan, downwash at the tail 
location is mostly given by the downwash caused by the wing vortices 
alone except near the junctiire. Figure 11(a), for an r/sm of 0.4, 
exhibits a decrease in downwash from the wing vortices, and the total 
downwash is not as well given by the wing vortices as in figure 6(a) 
especially near the body- tail juncture. Figure 8^), for an r/sm of 
0.6, shows that except near the body-tail juncture the downwash in this 
case is largely that given by the body upwash. Finally, figiire 11(b) 
for an r/sm of 0.72, showing the same trend as figure 11(a) and fig- 
ure 6(a), indicates that for very large ratios of body radius to wing 
semispan, such as for canard configurations, the downwash at the tail is 
almost entirely determined by the body upwash. 


Comparison of Rolling Up of Vortex Sheet Behind 
Wing Alone and Wing-Body Combination 


Plane wi ng and body .- It is interesting to compare the rolling -up 
process of the vortex sheet behind a wing alone with that behind a wing 
body combination under the conditions of sketch (c). 
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Both vings have the same aspect ratio, root chord, and Tq. However, 
this maximum T occurs at the center line for the wing (with elliptical 
distribution) and at the wing-body juncture for the wing-body combination 
(with near-elliptical circulation distribution) . Figure 5 of refer- 
ence 11 presents a picture of the rolling up of a vortex sheet behind an 
elliptically loaded wing as calculated by Westwater (ref. I6) in terms 
of the ratio of downstream distance d to rolling-up distance e. This 
ratio is related to p in figure 3 as follows (using eq. ( 5 ) of ref. 11 ): 




b wing span 

total lift coefficient 

K constant based on the shape of the trailing-edge circulation distri- 
bution 


or 


where 
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The value of K is 0.28 for elliptic loading. Now equation (lO) of 
reference 11 is: 


ksYo 


tcA 


where Pq corresponds to elliptic loading for the wing alone. When the 
last two equations are combined. 


d 

h 



( 0 . 28 ) 


4s 

^(^o Ao) 


or 


d ^ 2.24 f d^ s 
s ^ To/Vo 


In the present notation. 


d _ X 

s ^ 

The lift coefficient of a triangular -wing and infinite cylindrical -hody 
combination in terms of the wing alone lift coefficient is given by 
equation ( 38 ) of reference 15 as: 






It is clear that while the same Tq is chosen in this comparison, 
the total lift of the combination is less than that of the wing alone 
because of the loss of exposed wing area. For the calculations of fig- 
3) ^o/Yq equals 0.575 and Sm equals 3-75- These values being 
ass\xmed in the equation for d/s, 



For figure 3: 



0.2n 


Therefore the streamwise stations in figure 3 correspond to 0.2, 
0.4, 0.6, . . .2 wing semispans downstream. The values of d/e in 
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figure 5 of reference 11 correspond to values of x/sm of 0 , 0 . 233 , 
0-977^ 1*953^ 4.370, and 8.04o. Cross sections through the vortex sheet 
at these six stations have been drawn in figiire 12 as well as the nearest 
corresponding sheet configuration from figure 3* From figure 12 it can 
be seen that the rolling-up process of the vortex sheet behind a 
triangular-wing and small -body combination closely resembles that for a 
wing alone under the conditions of the preceding sketch. 

Banked wing and body .- The qualitative difference between the 
rolling-up motion of a vortex wake behind an isolated wing and wing-body 
combination both banked at a given small angle and at angle of attack 
can be reasoned simply with the aid of sketch (d), illustrating condi- 
tions immediately behind the trailing edge . 



Sketch (d) 

For the isolated wing, each panel vortex induces a velocity on 

the opposite vortex and both move in the dotted direction with uniform 
velocity. The motion in this case is symmetric with respect to the t)V 
plane. Now the addition of a body, with potential crossflow symmetrical 
respect to the qp plane, adds the same induced velocity at each 

vortex. The wing vortices in this case then move initially in the direc- 
tion of the resultant velocity Vp. Henceforth, the velocity due to the 
body crossflow is different at each vortex. Thus it can be seen that 
the wake motion is different from that behind an isolated wing, that is, 
completely nonsymmetrical. Similar reasoning for wing vortices of 
unequal, strength leads to the same conclusion. 

For a cruciform wing and body combination at 45° angle of bank and 
at angle of attack, the superposed velocity fields (from the wing and 
body) are both symmetrical with respect to the qp plane. Therefore 
the computed vortex wake motion should agree qualitatively with results 
for isolated cruciform wings. 
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II.- PRESENTATION AND EVALUATION OF METHOD 


The remainder of* this report will examine in detail the line-vortex 
method outlined in reference 13 . The manner by which the initially 
required circulation distribution was obtained is first discussed. The 
method of replacing the wake circulation distribution and the afterbody 
by discrete vortices then follows. The determination of the flow field 
in transverse planes and the stepwise integration technique are next 
explained. After a few sample calculations, an evaluation is made of 
the errors inherent in the stepwise integration and of the effect of the 
number of vortices used on the computed downwash distribution at chosen 
tail locations. 


Before the stepwise calculation of vortex paths behind the wings of 
wing-body combinations can be started, it is necessary to know the cir- 
culation distribution along the wing trailing edge. In general, the 
circulation distribution F(y) and the wing span loading are not equiva- 
lent. The conditions for equivalence are derived in Appendix A. The 
method of this report has been applied to triangular-wing configurations 
which have readily obtainable circulation distributions. However, the 
stepwise calculation method is applicable to configurations with any 
wing plan form. 

For a plane -wing and cylindrical -body combination, the span loading 
is equivalent to the circulation distribution and can be obtained by 
chordwise integration of the complete pressure coefficient (Ap/q)y. 
Equations for (Ap/q)y are given in the slender-body analysis of refer- 
ence 15. For such configurations, the circulation about the wing at any 
spanwise station y is given by linear theory as: 


The notation for this and the following equations is defined in fig- 
ure 13^ which shows the more general case of a noncylindrical body. 

Equation (ll) of reference 15 furnishes the wing loading coefficient 
for a plane-wing and cylindrical -body combination as: 


Circulation Distribution at Wing Trailing Edge 



( 1 ) 
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Substitution of equation (2) in equation (l) and integration yields: 

r(y) - - r 4 ) (3) 


When the body is not cylindrical, equation (2) is modified to: 



(U) 


The last term of equation (4) is due to the spanwise velocity and 
does not contribute to the circulation distribution. Hence one obtains 
the required circulation distribution by insertion of only the first 
term of equation (4) in equation (l) and integration. The result is 
identically equation (3)* This shows that the circulation distribution 
r(y) along the trailing edge of a plane wing on a noncylindrical body is 
dependent only on the cross section of the configuration at the trailing 
edge, as could be anticipated. The circulation is, in fact, the jump in 
the crossflow potential at the trailing edge and hence must agree with 
the slender-body result, equation (3)» 

For cruciform wings on a cylindrical body at angle of attack and 
sideslip, the complete velocity potential (eq. (47) of ref. 15) is: 


= 


■•¥{[■(>•$>■ 


^ cos 29 + ( 1 + 


r,->( 1 - I-:Y+ 


4r cos^ 20 + ( 1 + 


- 23= ( 1 +^ 4 ) (it cos ae 


1 1 
2 


VqU’z ± 


Vo3* 

*/2 


1 + ^ ] r^2 cos 20 + t2 ( 1 + ^ 


(1 . 4^ Y + 4r^ cos2 29 + + + 


2t2(^l + ^ 


1 + ) r,2 cos 20 


4 ; "1 


1,1 
2 1 ? 


+ Vop'y 


(5) 


NACA TN 3227 


17 


where 


and 




0 = tan”^ y* 


+ means 0 < 9 < n 
- means n < 6 < 2it 


The transformation to a configuration pitched a radians and banked 
cp radians is for small angles a: 

a' = a cos cp, p’ = a sin cp (6) 

Now for banked cruciform configurations, the span load distribution 
is not equivalent to the circulation distribution. Hence, instead of 
equation (l) one must use the general relation between F and the jtimp 
in the velocity potential at the wing trailing edge: 

r=Acp^ (7) 

For horizontal and vertical surfaces respectively, one obtains by 
substitution of equation (5) in equation (7): 

r(y) = - y^)(sm 2 y 2 - r^) (8) 

PY Q t 

r(z) = — 2 — ( 9 ) 

It is thus seen, by comparing equations (8) and ( 9 ) with ( 3 )^ that 
the circulation distribution for a cruciform wing and body combination 
is derivable by assuming that each wing acts independently of the other 
as part of a plane wing and body combination at an angle of attack given 
by equation (6). This can be seen also by the linear superposition of 
potentials in equation (5)* 

The foregoing equations are based on the slender-body theory of 
reference 15^ which is postulated for slender wing and body configura- 
tions at subsonic, transonic, and supersonic speeds. In order for the 
equations to apply at supersonic speeds, the entire configuration must 
lie well within the body nose Mach cone. Experimental data on models 
conforming to these restrictions agree well with the theory. In order 
to apply the results of slender-body theory to nonslender wing and body 
combinations, an aspect-ratio correction factor is employed. This factor, 
which for triangular-wing configurations is the ratio of the lift of a 
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triangular wing calculated by linear theory to the lift calculated by 
slender-body theory, is discussed fully in Appendix B. The results may 
be s\jimiiarized as follows: For subsonic -leading-edge wings (p tan w <l), 

equations (3)^ (8), and ( 9 ) should be multiplied by the factor 

k = i ^*/l - p^tan^ii) ^ 

For supersonic -leading-edge wings, (p tan to > l), these equations should 
be multiplied by equations (B?) or (BI 8 ) for the ranges of y indicated: 


k 



1 _ 

tan^ 


to 


Am - y 

V Sm + y 


(B7) 


for 


sm -r 

r + < y < Sm 

P tan to 
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1 

y(sm^ - y^)(3^ tan^ to - 1) 



+ 



(sm - y)sin"^ 


yP^ tan^ to - Sm 
P tan io(sm - y) 


(sm + y) sin"^ 


yp^ tan^ to + Sm 
ptan to(sm + y) - 


(BI 8 ) 


for 

r < y < r + 

P tan (0 


Replacement of Wake Circulation Distribution 
and Body by Discrete Vortices 


After the circulation distribution F(y) along the wing trailing 
edge has been determined, the next step is the replacement of the conse- 
quent wake circulation function by a finite number of vortices. A plot 
of the circulation function, equation (3), for example, looks like that 
portion of the solid cvirve above the wing in sketch (e): 
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The loading over the body, which can be calculated from the corre- 
sponding body pressure coefficient (eq. (12) of ref. 15), is not explic- 
itly required here. It will be seen that this loading is automatically 
accounted for by the vortex model set up for the wing loading. The 
curve r(y) over the wing can be replaced by a suitable number of step 
functions (dashed lines in sketch), each of which corresponds to a vor- 
tex of strength equal to the step height and located at the step abscissa 
as shown. Although an infinite number of step functions would be 
required to duplicate the given curve, in practice the number selected 
are the fewest which render possible a reasonably acciirate prediction of 
the trailing -vortex sheet behavior. Single-, three-, and ten-step func- 
tion distributions have been used herein. Quantitative results will be 
presented later concerning the number of vortices used to replace the 
circulation distribution. For a single-step (one vortex) approximation, 
integrating equation (3) for the area under the f(y) curve and equating 
the result to a single rectangular area lead to the simple result shown 
in figure l4, namely, that the spanwise distance from the body of a 
single vortex representing the circulation distribution on one panel of 
a plane, subsonic -leading-edge wing and body combination is approximately 
constant and equal to about O.76 times the exposed wing semispan for all 
ratios of body radius to wing semispan. In this approximation by a 
single vortex the value of at the wing-body juncture is assigned 

to the vortex. In all cases the two panels of the wing are treated inde- 
pendently. A circulation distribution which, unlike the preceding sketch, 
is asymmetrical across the wing trailing edge will, therefore, require 
an asymmetrical vortex distribution. 


The effect of a circular body in the presence of this wing vortex 
distribution is accounted for by placing an image vortex for each wing 
vortex within the body at the inverse point, as indicated by sketch (f ) . 
The image vortex cancels the velocity normal to the body due to the wing 
vortex . 
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Sketch (f) 


The foregoing discussion applies to the trailing -vortex filaments 
from the wing and the image vortices contained within the body. These 
filaments can be considered joined by bound vortices within the wing 
and adjacent body to form horseshoe vortices as shown in sketch (g). 

The portion of the total lift of 
the wing-body combination carried 
by the wing and by the body is rep- 
resented by the bound filaments of 
length Z2 within the wing and 
within the body, respectively. A 
plot of r(y) across the body is, 
therefore, obviated since the single 
vortex representing the body loading 
automatically appears at the inverse 
(image) point of the wing vortex. 
This vortex model, due to Lennertz 
(ref. 17), assumes that the ratio 
of body to wing lift is proportional 
to the ratio Z^^/Zg (Z^ + Zg being 
the "effective semispan") , and is 
known from experiments to represent 
the distribution of lift between 
wing and body with good accuracy, 
deviating but slightly from slender - 
body theory. 

Consistent with the assumed two-dimensional nature of the solution, 
the segments of the vortices bound within the wing and adjacent body are 
neglected. It is then assumed that the remaining filaments extend 
upstream to infinity. 

Determination of the Flow Field in Transverse Planes 



The nonlinear differential equations governing the downstream motion 
of the vortex wake require a stepwise solution for the vortex positions 
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at successive streamwise stations (or time increments). Consider a 
vortex pattern at any downstream station such as section A-A in figure 2. 
(The image vortices inside the body are not shown in the figure.) It is 
ass\omed that every vortex influences the flow field according to the 
Biot-Savart velocity law of planar, incompressible-flow vortex theory, 
as illustrated in the sketch (h). This assumption will be justified 
later . 



Vi = 


2irr5 


( 10 ) 


The vertical component, Wi, of Vi is: 

r, 


Wi e ^ 


(^i - Aj) + (’ll - 'Hj)' 
The horizontal component, Vi, of Vi is: 

-rj % - ^j 


Vi = rAi = 


(Ai - Aj)^ + (T)i - -)2 


( 11 ) 


( 12 ) 


In general, for a flov field containing n vortices, the velocity 
components induced at the position of the ith wing vortex by the other 
n-1 vortices are: 

^ Y ~ ^j^ 


( 13 ) 
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(lA) 


Equation (lO) is derived from consideration of an infinitely long, 
incompressible-flow vortex filament, and is, therefore, applicable only 
to two-dimensional flow. For supersonic speeds, the use of equation (lO) 
is consistent with the slender-body theory underlying this study for 
these reasons: Reference 7 shows that an infinitely long vortex filament 

parallel to a supersonic stream obeys the Biot-Savart law of incompres- 
sible flow. Omission of the bound vortex and wing chord loading causes 
the difference between initial and asymptotic downwash at the wing trail- 
ing edge seen in figure 15 reproduced from reference 7; the downwash 
along the wake center line approaches the asymptotic value within a chord 
length behind the trailing edge for the triangular wings of PA < 4. 

Hence the asymptotic downwash given by a two-dimensional trailing-vortex 
system has been used. 


Superimposed on the velocity field due to the wing and image vortex 
system is the body potential crossflow velocity field. From the poten- 
tial: 

cPi =V^ ^ ) ( 15 ) 

Vy^ + z* y 


where z' is measured from the body axis (instead of the wind axis), one 
readily obtains the velocity components induced by the body: 
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or in dimensionless form: 


w 



(A^ + 


^ - -2a An' 

^O (A^ + 


( 16 ) 

(17) 


(18) 


(19) 


NACA TN 3227 


23 


where 


T]' = T] + ^ tan a * Ti + |ia (20) 

This transformation from the usual potential equations is neces- 
sitated here hy the fact that the coordinates are not "body centered, and 
the hody is inclined, with respect to the free-stream direction. 

These last equations assume that the body angle of attack is small 
enough so that tan a ~ a, cos a 1 (which is true for a up to approxi- 
mately 20°), and that the contribution of the crossflow to the local 
streamwise velocity is negligible. As stated previously, viscous cross- 
flow separation around the body, with its consequent vortex wake, is not 
considered here. 


Stepwise Determination of Vortex Paths 


The replacement of the wing trailing-edge circulation distribution 
by step functions and associated vortices, together with the placement 
of image vortices within the body, as detailed above, provide the start- 
ing point for the stepwise calculations to be described. The spanwise 
and heightwise coordinates of all vortices at p = 0 are known. The 
downstream incremental motion of these vortices will next be determined. 

Now at any point in a transverse (yz) plane, the lateral velocity, v, 
is given by the sum of equations (ll<.) and (19). The vertical velocity, w, 
is similarly given by the sum of • equations (13) and (18K The streamwise 
velocity, u, perpendicular to the plane, is everywhere Vq. Consider 
the point where a vortex line passes through the transverse plane. In 
a time interval At the vortex filament moves with the fluid a vertical 
distance Az = w At. Since At = Ax/Vq, Az = (w/Vq)Ax, then nondimen- 
sionally: 

Aq = ^ Ap (21) 

''o 

Similarly the dimensionless lateral movement is: 

AA = ^ Ap ( 22 ) 

''o 

The last two equations then furnish points in transverse planes a 
distance A^ apart, the locus of which is the path of a vortex filament 
moving in the wing-body flow field. With each vortex, a "strength" F^/v 
(rather than fj) will be associated in order that equations (13) ^ ° 

and (Ih) will actually yield the velocity ratios w/Vq and v/Vq. 
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In the plane p = 0, which contains the wing trailing edge, as long 
as the wings have no incidence with respect to the body, all the vortex 
filaments have coordinates r| = 0. Therefore, the total spanwise veloc- 
ity induced at the location of any given wing vortex can readily be seen 
to be zero, from equations (l^i-) and (19)* In fact, the initial motion 
of such a vortex sheet is always one of pure vertical displacement. The 
total vertical velocity summed from equations (13) and (18) is multiplied 
by a suitable increment Ap yielding an incremental distance Aq moved 
by the given vortex line, according to equation (2l) . Values of Aq are 
thus calculated for all the wing vortices, the entire group of which is 
so transported to the plane p = Ap. The coordinates of each wing vortex 
in this new plane are found simply by adding the computed increments to 
the coordinates at the previous station. Of course, symmetrical proper- 
ties are used wherever possible to obviate the calculation of the paths 
of vortices from each half of a wing. Corresponding to the new location 
of each wing vortex line at the station p = Ap, the image vortices 
within the body are repositioned according to the following formulas: 


A 


image 


+ (q + p tan a)‘ 


(23) 


q + p tan a 

q. „„ = -p tan a + 

'image ->2 / j. \2 

A + (q + p tan a) 


(24) 


where (A^t]) are the coordinates of the wing vortex at the new station 
]i = A)-l. This readj-ustment of the image vortex positions, required hy 
the downstream displacement of the inclined hody away from the free- 
stream direction, is made because the image vortices are hound within 
the hody; the image vortex displacement, in fact, measures the amount 
of lift carried hy the afterbody (see ref. 13)* 


With the new positions of wing and image vortices this procedure is 
repeated. That is, the total lateral and vertical velocity induced at 
each wing vortex location is calculated and multiplied hy an incremental 
distance A[i. The new wing vortex positions are obtained hy adding the 
computed increments to the original coordinates. The image vortices are 
then repositioned according to the new wing vortex locations, and the 
stepwise calculation continued to the desired do^mstream station. 

This stepwise procedure is simply a method for integrating numeri- 
cally ' n simultaneous differential equations of motion of n/2 wing 
vortices for the downstream paths. A closed solution for the three- 
dimensional paths is, in general, extremely complicated, if not impos- 
sible. While the equations of the vortex motions in a transverse or 
crossflow plane can he written, introduction of time dependence (equiva- 
lent to the streamwise coordinate) renders a closed solution extremely 
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difficult. Reference 13^ using complex variable notation, cites explic-- 
itly the entire path of a pair of vortices in the presence of a cylinder 
in the transverse plane. The equations vere derived hy Villat in refer- 
ence 18 . Even for this simple case, where n = U, recourse is necessary 
to numerical solution for the time or streamwise coordinate. This solu- 
tion will he analyzed in greater detail subsequently, as well as the 
manner of choosing the proper spacing Ah of successive stations. 


Sample Calculations 


Three examples will now be presented to illustrate the calculative 
procedure and resulting vortex paths. They are: 

1. A plane -wing and cylindrical -body combination at angle of 

attack. 

2. A cruciform- wing and cylindrical-body configuration at angle 

of attack and bank. 

3 . A plane-wing and cylindrical-body combination with the body at 

‘angle of attack and the wing panels at differential incidence. 

The procedure has been systematized by prepared computation forms 
such as table III which is used when the circulation distribution on one 
wing panel is replaced by a single vortex. 

Example 1, plane triangular -wing and cylindrical -body combination 
at angle of attack . - The data for this example are: wing aspect 

ratio = 2/3, ratio of body radius to wing semispan = 0.6, angle of 
attack = and Mach number = 2.0. 

The leading edge of the wing is subsonic (3 tan w = 0.2912 < l), so 
the trailing -edge circulation distribution r(y) is calculated from equa- 
tion ( 3 ) and multiplied by the aspect-ratio correction factor k = 0.917 
(see Appendix B) . A plot of r(y) vs. y like that of figure lU is 
obtained, with a maximimi ordinate R/Vq = 0.12796 in this example. 

First, this distribution is replaced by a single vortex per wing panel, 
with the location consequently given by figure lU as s'/(sui-r) = 0.763^ 
or = 1.509 (r/sm = 0.6). The vertical coordinate of the vortex at 
the wing trailing edge (p = O) is zero in this case. The dimensionless 
spanwise coordinate of the image vortex is, therefore, = lAj “ 0.663» 

Referring now to table Ill(a), the above values of ^ 

Ag, and filled into the proper boxes of columns (^, (^, and By 

symmetry, the values of A^, A^, q^, and q^ can be readily filled in, and 
an increment Ap is chosen (l here). Columns ^ through based on 
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equations (13) and (l^)^ are then computed in order, the last two giving 
the vertical and horizontal veloci^ components induced at vortex 1 by 
the other three vortices (columns @ and respectively). Rows 
through based on equations (l8; through (22), next add the potential 
crossflow velocity components at vortex 1 . Thus, rows and give 
the new coordinates of vortex 1 for use in table Ill(b), while uhe 
bottom calculation furnishes the coordinates of the corresponding new 
image vortex 2 position. 

Columns through (J) of the table Ill(b), for station p = 1 , can 
now be filled in and the procedure repeated, as illustrated. This step- 
wise computation is continued for as many stations downstream as desired, 
each part of this example representing positions at increments of one 
body radius (Ap = l) downstream of the wing trailing edge. The results 
are sunmiarized graphically by figure l6. 

For the same given data, the calculated circulation distribution 
for one wing panel is now replaced by three vortices of equal strength 
F/Vq = O.OU2653. Their spanwise locations are determined graphically, 
equating areas under the curve of r(y) vs. y, and in this example are: 

= 1.652, Ag = 1.558? a-nd A3 = I.316. The three image vortex loca- 
tions are again calculated as the inverse points, giving: 'K = 0. 75988, 

As = O.6U185, and Ag = 0 . 60533 - The results of this case are shown in 
figure 17. 

Example 2 , cruciform triangular-wing and cylindrical-body configura- 
tion at angle of attack and bank .- The configuration data of the pre- 
ceding example are again assumed in addition to which another pair of 
wing panels now exists, forming a cruciform wing-body arrangement banked 
45°. The strengths of the assumed four vortices originating from the 
four wing panels are identical, as required by equations (8) and (9) • 
These equations with equation (6) provide the required circulation dis- 
tribution r(y) and f(z) which is corrected for aspect ratio as explained 
in Appendix B. Figure l 4 again furnishes the initial positions of the 
four wing vortices, whose corresponding four image vortices are then 
located at the inverse points. The calculative procedure follows that 
described above. Because of planar (xz) symmetry, only the paths of a 
pair of vortices on one side of the body need be computed. The results 
of this case are presented in figure 18, which shows that there is a 
tendency toward "leapfrogging," sketch (i )3 although the process appears 
to be retarded initially by the body. 


®The phenomenon of "leapfrogging" vortices is discussed in detail 
for isolated wings in reference 12. 
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Pair of vortices from 
upper panels 


\ / 
\ / 


\ / 
\ / 



Pair of vortices from 
lower panels 


ft 


Leapfrogging" paths in crossflow plane 


Sketch (i) 


Example 3 } plane -wing and cylindrical -body combination with the 
hody at angle of attack and the wing panels at differential incidence . - 
The given data for this case are the same as for the first example with 
two exceptions. One wing panel is deflected and the other is at zero 
incidence with respect to the hody axis such that the vortex from each 
panel has the same strength. Now the vortices from both wing panels 
are rotating in the same direction. The vortex shed from the deflected 
panel is assumed to start above the xy plane because of the angular 
displacement about the hinge line of that panel's trailing edge. It is 
further assumed that the vortex from each panel is shed at the same span- 
wise station. 

Figure I9 illustrates an effect characteristic of wings differen- 
tially inclined. The vortex from the wing panel which has no incidence 
to the body moves essentially the same way as in the plane wing case. 
However^ the vortex from the wing panel which is at a negative angle of 
attack travels sharply upward and inboard. This marked motion is due 
to two effects: First, the vortex starts at the trailing edge with some 

Ti > 0 because of the negative angle of attack of the panel, and second, 
the body crossflow and image vortex- induced velocities are ciomulative in 
this case instead of being subtractive as in the plane case. At a suf- 
ficiently high angle of attack, it may be possible for the rising vortex 
to jump over the body to the same side as the other vortex. Such an 
occurrence would cause a sharp increase in rolling moment from the tail 
surfaces. 

For the case of a pair of wing panels inclined together, the trail- 
ing vortex paths can be estimated qualitatively by consideration of the 
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vortex picture in the crossflow plane. Thus^ for a pair of wing panels 
inclined at a positive angle of attack on a "body whose axis is parallel 
to the free stream^ the vortex paths trail below the free-stream direc- 
tion with concavity upwards since there is no body upwash and the effect 
of the image vortices is to depress the wing vortex wake. For a pair of 
wing panels at a negative angle of attack on a body at a positive angle 
of attack^ the trailing vortices are inclined upward since the body 
crossflow and vortex- induced velocities are additive. 

These and all other computations except those with 10 vortices per 
wing panel were performed with a desk calculator. It was found that the 
solution of vortex paths downstream of the wing could be accomplished in 
about fifteen minutes per station for a single -wing vortex^ and two hours 
per station for a three-wing vortex scheme. In general, the number of 
computations per solution increases approximately as the square of the 
number of vortices assigned per wing panel. 


Evaluation of the Effect of Stepwise Integration 
on Vortex Paths 


In the present line-vortex method, the principal factors governing 
the labor expended in the solution to a given problem are the size of 
the "time" increment and the number of vortices used to replace the 

trailing vorticity. The first factor will now be discussed in detail. 

The necessity for a stepwise solution for vortex paths stems from 
the difficulty of integrating n simultaneous differential equations of 
vortex motion (eqs. (I3) and (l4), together with eqs. (18) and (19))* 

The differential equations are therefore solved with a small, finite Am*. 
In this report, the use of equations (2l) and (22) implies the use of 
Euler’s linear integration method, that is, 

(''i) A 


where ^ 

Vi = or T)j^, and i = 1, 2, 3^ • • • ^72 
A Taylor series expansion shows the higher order terms neglected: 
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Thus, it may be seen that the use of Euler's method without appre- 
ciable error requires either the second and higher derivatives to be 
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small or A|i to te small. The following discussion indicates the impor- 
tance of the higher order terms for vortex motion. 

Comparison of known exact solution of vortex motion with stepwise 
j.ntegration . - A simple example of the mutual interaction of two vortices, 
for which the exact solution is known, furnishes an insight into the 
nature of the cumulative error incurred "by this stepwise -integration 
method. It is easily shown hy the methods of hydrodynamics (e.g., p. 320 
of ref. 19) that two vortices of strengths and fg will, due to their 
own mutual influence, each rotate about their common ^centroid^^ in a cir- 
cular path with constant angular velocity which is equal to ( + F, ) /2jt 
divided by the square of the distance between the vortices. This motion 
is illustrated in figure 20. The solid lines radiating from the centroid 
represent the ends of a constant time increment (proportional to AM-) and 
the solid, spiral-like curves are the corresponding paths computed by a 
stepwise approximation. The dashed radial lines and spirals relate to a 
time interval half of that used before. It appears from figure 20 that 
the error or discrepancy between circular and spiral paths increases 
while the rate of growth of this error decreases with time (or distance 
downstream of the initial position). Further, it is seen that a decreased 
size of time increment results in a decreased error, for this example, in 
approximately the same ratio. Of course, this case of completely circular 
vortex motion cannot be solved accurately by Euler’s linear scheme without 
using extremely small increments. However, for many of the computed vor- 
tex paths behind the configurations considered in this report it will now 
be shown that Euler’s method is sufficiently accurate. 

Effect of size of increment Am on computed vortex paths . - The 
error inciirred by using the linear integration method can be judged by 
calculating vortex paths with different increments AM and extrapolating 
the results to AM = 0. This has been done for single- and three -vortex 
schemes at a = 5^ and typical results appear in figures 21 and 22. 

Both figures show that the stepwise method employing finite increments AM 
yields results which converge fairly rapidly to the exact solution (for 
which Am >dM — >0). Here the error, using a given value of AM^ is 
indicated by the difference between the value of A or t] extrapolated 
to Am = 0 and the value of A or t) at the given value of Am* As in 
figure 20, the error is seen to increase with time or distance m down- 
stream of the wing. For a single-vortex scheme (fig. 2l), the error is 
approximately proportional to pAM and is small, in general. The error 
becomes greater and less predictable, for the same increment size, with 
a multiple -vortex scheme (fig. 22). A single-vortex scheme is therefore 
less subject to integration error than a multiple -vortex scheme. The 
reason why the error increases nonlinear ly for the three -vortex scheme 
of figure 22 can be explained by reference to figure 22(c) which shows 
the projection of the vortex paths on the crossflow plane. For small 
values of p the vortex paths are fairly linear, with the vortex posi- 
tion error given by figure 22(a) being small and essentially linear 
with AM. Further downstream, as the effect of the rolling-up process 
becomes marked, the paths are seen to become spiral-like. Each vortex- 
position error curve becomes increasingly nonlinear as the vortex path 
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approaches a maximum or minimum vertical position in the crossflow plane. 
The use of smaller increments A^i then reveals the tendency toward con- 
vergence of the path to one of a lesser radius of curvature, as illus- 
trated previously for the case of truly circular paths. 

From these results it appears that the selection of a suitable 
increment A^l for any given case should he made hy trial calculation 
with several values of A^i for a few stations downstream of the wing. 

As long as the computed vortex paths are relatively linear in the cross- 
flow plane, the largest value of should he chosen which permits the 

extrapolated error in vortex location to he the maximum tolerable at the 
furthest downstream station at which downwash is desired; Ai-i should he 
decreased locally wherever the vortex paths appear to he approaching a 
maximum or minimum height in the crossflow plane. 

More accurate numerical integration schemes than Euler’s method can 
he found in reference 20. Curve -fitting and extrapolation formulas are 
presented which hy taking account of the higher order derivatives of the 
path enable one to maintain sufficient accuracy of vortex positions with 
a given A^I even when the paths are markedly nonlinear. 

Comparison of computed paths with exact solution for vortex paths 
in crossflow plane .- The exact solution known for the vortex path in the 
crossflow plane corresponding to the single -panel vortex representation 
mentioned previously can he used to check the accuracy of the stepwise 
integration. Equation (VI-4) of reference 13 cites the result obtained 
hy Villat (ref. I8) . In the present notation, it is: 
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At an infinite distance downstream the vortices are infinitely far 
removed from the body (t] + pa — >00) so that the asymptotic spacing 
2(A)^_^ of the vortex pair is given hy the relation: 

2(^)u=« - In 2(A) = 2 Aq f 1- In 
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Two solutions by the stepwise method have been checked against 
equations (25) and (26) • This has been done by inserting pairs of values 
of the coordinates q. and t), from the stepwise -calculated paths^ into 
equation (25) and solving for the A coordinate by trial. The compari- 
son appears in figure 23 , showing the results to be practically coincident. 


Effect of Number of Vortices on Downwash at Tail Locations 


Comparisons were, made of the downwash along a line, representing a 
horizontal tail plane, resulting from the separate presence of one, 
three, and ten vortices, all derived from a given circulation distribu- 
tion. Possible tail locations for both a missile-type and high-tail 
airplane -type configuration were chosen. 

Missile-type configuration .- For the former, |i = 10 was selected; 
that is, a tail location 10 body radii downstream of the wing trailing 
edge. Therefore downwash was first computed along the line p = 10, 
q = - 10 tan a, 1 < A < 6, which is a horizontal line through the body 
center line of the missile -type configuration. 

Figures 24(a) and (b) show the effect of number of vortices on down- 
wash at the chosen missile tail location with and without components of 
the body potential crossflow in the stepwise computations. These omis- 
sions were made for two reasons: One was to determine the extent 'to 

which -the stepwise calculations could be simplified -without obtaining 
dissimilar results in do'wnwash, and the other ■was an attempt to account 
for the flow around the body at high angles of attack. As an assumption 
for> the high -angle -of -attack crossflow, the potential crossflow (eqs. (18) 
and (19)) was replaced by equation (I8) alone with q* = 0. This assump- 
tion implies that the velocity distribution on the lee side of the body 
is the same at any vertical station as along the horizontal diameter, as 
indicated in sketch (j). 



Potential crossflow. 
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For the aspect ratio 2/3 configuration of figure 2^i(a), with 
r/sm = 0.6, the numher of vortices used has a negligibly small effect on 
downwash except when the body crossflow is entirely omitted in the step- 
wise calculations. Observance of the various vortex positions at the 
tail location, shown at the top of figure 24(a), shows that the vortices 
from the wing panels are subject to a relatively strong body crossflow 
so that their computed positions are greatly in error if the body upwash 
is neglected. Thus, when at least the upwash component of the body cross - 
flow (the center line crossflow in sketch (j)) is added to the stepwise 
calculations there is no appreciable effect of either the number of vor- 
tices used or the strength distribution for a given multiple -vortex 
scheme. Figure 24(a) also shows that only a small change in downwash 
results from replacing the potential body crossflow by the assumed center 
line crossflow. 

For the aspect ratio 4 wing and body combination of figure 24(b), 
with r/sjn = 0.2, it is seen that there is little effect of body cross- 
flow on the downwash at the tail. This could be expected because of 
the relatively large wing. While the use of a single vortex does not 
adequately approximate the magnitude of downwash over a tailspan equal 
to the \>ringspan, a 3“Vortex scheme appears to be as satisfactory as a 
10-vortex scheme. Reference to the computed vortex locations shown at 
the top of the figure indicates that the vortex sheet from a high-aspect- 
ratio wing at low angles of attack rolls up so slowly and extends over 
such a large spanwise distance that a single vortex cannot give the cor- 
rect distribution of downwash. 

At higher angles of attack, the vortex wake will be even further 
from the horizontal tail plane because the body will be inclined below 
the free -stream direction more than the vortex wake according to fig- 
ure 7* Although the vortex strengths increase directly with angle of 
attack, the downwash at the tail decreases very rapidly with distance 
from the wake. Consequently, the effect of number of vortices should be 
smaller at higher angles of attack. 

Airplane -type configuration . - Now a possible horizontal tail loca- 
tion for a high-tail airplane-type configuration is p = 3 and t) 2 2, 
that is, 3 body radii doimstream of the trailing edge of the wing and 
2 body radii above the body axis. Figure 25 presents the effect of the 
number of vortices used on the downwash at this tail location. The 
aspect ratio 2/3 configuration of figure 25(a), with r/sm = 0.6, shows 
no noticeable effect of number of vortices on either the distribution or 
magnitude of downwash simply because of the distance of the horizontal 
tail plane above the vortex wake. Neither is there any appreciable effect 
of assuming the body crossflow to be the center line crossflow. The 
aspect ratio 4 configuration of figure 25(b), with r/sm = 0.2, shows 
that 3 vortices give the same downwash as 10 vortices. However, from 
figure 7 the vortex wake behind the wings of both high- and low-aspect- 
ratio conf igiirations at small values of p is seen to be fairly close 
to the free -stream direction at all angles of attack. Therefore it can 
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be expected that with increasing angle of attack, the tail plane 
approaches the wing vortex wake and the use of at least three vortices 
to represent the wake behind wings of such high-tail airplane-type con- 
figurations will be required to obtain the correct magnitude of downwash 
across the tailspan. 

To obtain the shape of a vortex wake behind a wing, clearly a mul- 
tiple trailing-vortex system, such as the 10-vortex systems used herein, 
is required. However, for the computation of downwash at tail locations 
a minimum number of vortices is desirable and this number depends mainly 
on the distance between the wing wake and tail surfaces. The number 
varied from about one to three for the configurations studied herein, 
depending on whether that distance was large or small. An empirical 
relation, based on the calculations made in this study, for estimating 
the nearest integral number of vortices required per wing panel is: 


where h is the height of the horizontal tail above the body axis in 
terms of body radii. 

Note on calculations involving large numbers of vortices .- In set- 
■ting up the initial (trailing-edge) vortex distribution from the given 
circulation distribution r(y) , one must exercise great care when n is 
large (i.e., n > 20) . The spacing of the vortices which approximate 
r(y) stepwise must be checked to insure a smooth curve of their divided 
differences. Otherwise, spurious effects such as the "loss" of a vortex, 
caused by excessive induced velocity, sketched below, can occur (see 
fig. 24(a)). 


N = 1 + 


0.16a 




Tf 


Lost" vortex 



Initial vortex distribution 
at n = 0 


Incorrect result at n = A|i 
due to vortex No. 2 


Sketch (k) 


CONCLUSIONS 


An investigation has "been made of a line -vortex method for computing 
vortex paths^ downwash^ and sidewash behind wings of wing -body combina- 
tions. Application of the method to the calculation of wake shapes, 
vortex paths, and downwash behind triangular -wing and cylindrical-body 
combinations with various prescribed data leads to the following con- 
clusions; 


1. Wake shapes calculated by two-dimensional line -vortex theory 
agree qualitatively at least with wakes observed in a supersonic wind 
tunnel. Calculated vortex paths using this line-vortex method agree 
well with a known exact crossflow-plane solution. 

2. A missile -type wing -body-tail combination with the horizontal 
tail located in the body diametral plane 10 body radii downstream of the 
wing trailing edge at 5^ angle of attack was considered. It was found 
that if the plane wing is of low aspect ratio (order of 2 / 3 ), downwash 
can be computed at the horizontal tail from the wing vortices alone for 
small values of the ratio of body radius to wing semispan (order of 0.2) 
and from the body upwash alone for large values of that ratio (order 

of 0.7)- For high-aspect-ratio wings (order of 4) on small bodies, do\m- 
wash at the tail location can be well approximated only by considering 
all the flow components - wing vortices, image vortices, and body cross- 
flow. 


3 . A multiple trailing-vortex system is, of course, required to. 
determine the shape of the wake behind a wing. However, for the computa- 
tions of downwash at a tail location the number of vortices required per 
panel depends mainly on the distance between the wing wake and tail sur- 
faces. The number varies from about one to three for the configurations 
studied herein, depending on whether this distance is large or small. 

An empirical relation, based on the calculations made in this report, 
for estimating the nearest integral number of vortices N required per 
wing panel is: 


N = 1 + 


0 . 16 a 


r 



h 


- \ia 


where h is the height of the horizontal tail above the body axis, \i 
the distance of the tail from the wing trailing edge (both in terms of 
body radius), a is the angle of attack in radians, A is the wing aspect 
ratio, and r/sm the ratio of body radius to wing semispan. 

4. The rolling up of a vortex wake behind wings of unbanked wing- 
body combinations qualitatively resembles the wake pattern behind com- 
parable isolated wings. The comma-shaped rolling-up pattern of the 
vortex wake behind a subsonic-leading-edge triangular wing panel, with 
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the eventual single-vortex core, is due to the nearly elliptic circula- 
tion distribution along the wing trailing edge. In contrast, a supersonic- 
leading-edge triangular wing panel generates an S-shaped vortex wake 
behind the portion of* the trailing edge which has a linear circulation 
distribution, and the S-shaped wake eventually rolls up into two vortices. 

5. The spanwise distance from the body of a single vortex repre- 
senting the circulation distribution on one panel of a plane, subsonic- 
leading-edge wing and body combination is approximately constant and 
equal to about O.76 times the exposed wing semispan for all ratios of 
body radius to wing semispan. 

6 . The type of vortex wake, and resultant downwash from a wing, is 
significantly affected by the circulation distribution along the wing 
trailing edge. It should be noted that the circulation distribution 
which must be used in setting up the initial vortex distribution along 
the wing trailing edge is not equivalent to span loading when the wings 
are in sideslip or the body noncylindrical. 


Ames Aeronautical Laboratory 

National Advisory Committee 
Moffett Field, Calif., 


for Aeronautics 
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APPENDIX A 

NONEQUIVALENCE OF WING SPAN LOADING AND CIRCULATION DISTRIBUTION 
r(y) ALONG THE TRAILING EDGE OF A WING 

The circulation function F in fluid dynamics is defined by the 
contour integral: 


Let r be evaluated around a streamwise chord of a flat plate such 
as AB in figure 13 . Then: 


From the Bernoulli equation one obtains the lifting pressure coef- 
ficient for a general wing -body combination as: 



where 


q = ui -I- vj + wk 


and 


ds s dxi + dyj + dzk 


r = 


I 


u dXj^ 



o 

u^ dxj^ + u^ dxi 



where u^i means (u)^. 


where means (u)^ >+o^ means (u)^ >-o* 

Within the limits of linear theory^ u^ = -u^, so that: 



or 



(Al) 


The span loading on a wing is defined by the relation: 



(A2) 



(A3) 
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For a noncylindrical body-wing combination, two linear crossflow 
velocity potentials must be added; associated with a cylindrical 

body-wing combination and expanding body can be superposed 

as follows: 




When the potential subscript 2 is dropped. 


^2 - 'iu = f('P) = ■^2 " ■^u = S(^) = = 0 

^2 = l>2l I 'lx I 

^2^-^u^ = 2lv<j^|(|V(pJ + |v^|) = 1^1 II vcp^l 

using the symmetrical properties of each potential. 
Substituting equations (a 4) into (A3), one has: 


(a4) 


_ 4ucp 


4vj)Vcp 


(A5) 


where the derivatives are evaluated on the wing (z=0). Equation (4) in 
the text is the result of substituting the expressions for potentials 
cp and $ in equation (A5) • It can thus be seen that the second term in 
equation (4) representing the contribution of the spanwise velocities 
(the second term of eq. (A5)), can be of the same order of magnitude as 
the first term, and not negligible as in the linearized theory for wings 
alone . ^ 


^The above derivation appears in reference 21 where it is pointed 
out that for axisymmetric bodies, crossflow velocities may considerably 
exceed axial velocities, that is, 0(u) < 0(v). 
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Now the span loading is obtained by substituting equation (A 5 ) in 

(A2): 


8y 



dxj^ 



dx 


1 


and using equation (Al) : 






(A6) 


Thus it is seen that the span loading along the wing trailing edge 
of a noncylindrical body-wing combination differs from the circulation 
distribution by the last integral term. In fact, it can be shown, by 
writing equation (A3) in terms of velocity potential and comparison with 
equation (7), that to any order theory span loading is not equivalent to 
circulation distribution as long as sideslip angle of flow is not zero 
or the body noncylindrical. For a banked cruciform configuration, there- 
fore, this nonequivalence is always present. 
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APPENDIX B 

CORRECTION OF SLENDER -BODY THEORY TO ACCOUNT 
FOR MODERATE ASPECT RATIO 


The theory developed in reference 15, for the lift of wings and 
wing -body combinations, is based on the assumption of slenderness. The 
aspect ratio of the configuration analyzed by this theory must be low 
enough so that the wings are near the axis of the Mach cone. A correc- 
tion factor is required, therefore, to apply slender-body results to 
higher-aspect-ratio wings according to the formula: 


S/+B ~ ^^^^^+B^sbT 


(Bl) 


The correction factor k is obtained by comparison of slender-body 
theory with the more exact linear theory (applicable to high-aspect-ratio 
wings) for triangular wings (ref. 22), that is: 


k = 


(Ly). 


LT 


(Bw), 


SET 


(B2) 


It has here been assiamed that this wing aspect-ratio factor k can 
be used for the winged portion of wing-body combinations and for other 
than triangular -wing plan forms. Thus one can correct the lift of a 
higher-aspect-ratio configuration by means of equation (Bl) . This assump- 
tion is justified for small values of the parameter |3 tan w (see fig. 26) 
when k is approximately one, and slender-body theory is itself applica- 
ble. Satisfactory results should also be obtained for large wing-to-body 
area ratios since the wing then carries the major portion of the lift of 
the combination. Experimental data confirm the validity of these assump- 
tions for wing-body configurations of aspect ratios on the order of 
1 to B. 


For triangular wings with subsonic leading edges, the lift distribu- 
tion given by linear theory has the same shape as that given by slender- 
body theory. Therefore the k factor can be defined as in equation (B2) 
by total lift ratios, and is independent of the spanwise coordinate y. 

For such wings, the method of reference 15 (applied to wings alone) yields 
results identical to the low-aspect-ratio triangular-wing results of 
Jones (ref. 23 ). However, it is known that Jones' results overestimate 
lift-curve slope when the parameter 3 tan to is not small relative to 
unity. The factor by which lift-curve slopes calculated by slender-body 
theory must be multiplied to agree with linear theory is given as: 


4o 
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k = ■■ ■■■ ■ p tan CO < 1 

eL/ 1-P^ tan^ CO j 

, 2 

k = P tan (0 > 1 

Jt p tan (0 

These equations have been plotted in figure 26. 


1 


J 


(B3) 


However, for supersonic leading edges (p tan co > l), the lift dis- 
tributions are not functionally identical so that equation (Bl) must be 
written in terms of local chord lift. Then k becomes a function of y, 
namely, the ratio of section lifts or circulation by linear and slender- 
body theoiy . Now in linear theory, the pressure coefficient on a 
supersonic -leading-edge triangular wing outboard of the Mach line is: 


q 


4a 



tan^ A 


(b4) 


Substitution of this expression in equation (l) yields the circu- 
lation distribution (in this case equivalent to the span loading): 


r(y) = 


tan A(sm - y) 
^ - tan^A 


(B5) 


For a wing alone r = 0 in equation (3), thus giving the slender- 
body theory result: 


r(y) = - y^ 


(B6) 


Hence, for the part of the wing trailing edge outboard of the Mach 
cone, the factor k for a wing -body combination is the ratio of equa- 
tion (B5) to (B6), or: 


k = 


sm - y 


J P^ tan^ (0-1 J Sjn + y 


(BT) 


for 


Sm - r 

r + < y < Sm 


P tan to 
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Fortunately, the accuracy of the k factor can he checked in this- 
case since an exact answer for r(y) for the wing -body combination can 
be calculated by linear theory without excessive labor. Referring to 
figure 27, at any point P(x,y) on the wing trailing edge outboard of 
the point of intersection of the Mach cone with the trailing edge, the 
potential cpg can be evaluated by integration of the distribution of 
sources at points (a,r) within the forecone of P, that is: 





JT 


If 

forecone 
Of P 


(B 8 ) 


where, in the plane of the wing (z = O), the hyperbolic radius r^ is: 


rn = y (x - a)^ - | 3 ^(y - t)^ 


(B9) 


and the angle of attack a is the sum of the body angle of attack, 
and the body angle of upwash: 


a = a-Q\l + — 


(BIO) 


On the wing, the forecone of point P produces the shaded area 
shown with the corresponding boundaries of integration. Thus: 


9g(x,y) = 




dT 


3 y-x 


x- 3 (y-r) 


1 + 72 )do 


T tan A y (x-a)^ - 3^(y-T)^ 


3 -tan A 


3 y+x 
3 + tan A 


^x+p(y-T) 


dT 




1 + ) da 




T tan A y (x-a)^ - 3^(y-T)^ 


(Bll) 


The integrals not containing r have been evaluated in refer- 
ence 24 , and in the present coordinate system give: 


y 3 ^ - tan^ A 


(B 12 ) 
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The integrals containing r^/rS yield: 


^5 = 


VoCCBr^(x-y tan A) 


y^/Py^ - 


(B 13 ) 


Therefore the total potential at P is 9 ^ = cp^ + cp^ or: 


93(x^y) = - VQag(x-y tan A) ^ 


1 r" 


J ^ - tan^ A y\/ P^y^ - 


(Bl4) 


The circulation distribution r(y) along the wing trailing edge 
X = cr is then (see eq. (?)): 


r(y)Tg = 2Voag tan ACs^-y) 


y - tan^ A y\/ 3^y^ - Sj^^ tan^ 


A 


where 


(BI 5 ) 


Sm - I* 

r + — - — tan A < y < Sm 


Thus equation (B15) furnishes the linear-theory answer for r(y) 

along the trailing edge outside of the Mach cone, while equation ( 3 ) 

multiplied by equation (BT) is the approximation to the linear theory 
assiamed by equation (Bl) with k = k(y) . This approximation is com- 
pared with linear theory in figure 28 at the spanwise position of the 
Mach line y = r + (sm-r)/ptan to for various values of r/sm as a 

function of 3 tan to. It is seen that the approximation asstimed by 

equation (Bl) with k = k(y) is satisfactory, producing a maximum error 
of about 20 percent and only about 10 percent for the range of the 
parameter 3 tan to from 0 to 4. 

For the region of the wing inboard of the Mach cone, the pressure 
coefficient given by linear theory is (for a wing alone): 


q 


4a 


y 


- tan^ A 


1 I 3 ^ cot A - 1 

1 + i sin"^ 

7T 


3 ( cot A - ^ 


— 3^ cot A + 1 
i sin- ii, ^ 

^ 3 Aot A + 


(BI 6 ) 
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Integration of equation (l) then results in: 


r(y),pg 




tan^ to - 1 




(sm - y)sin"i 


yP^ tan^ to - sm 
( Sm " y) P tan to 


(sm + y)sin"^ 


yP^ tan^ to + Sm 3 


(sm + y)P tan to 


]} 


(BIT) 


■nierefore the factor k(y) to he used for r(y) along the wing 
trailing edge inboard of the Mach cone is the ratio of equation (BI7) to 
equation (b 6) or: 


k(y) 


_1 

a/ - y^)(P^ tan^ to - 1) 



(sm - y)sin 


yp^ tan^ w - 
(sm -y)P tan to 


(sm + y)sin"i 


yP^ tan^ to + sm 


(sm + y)P tan 


to -> 


(B18) 


for 


Sm " ^ 

r < y < r + 

P tan to 
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TABLE I.- VORTEX COORDINATES 


A, T] VERSUS DOWNSTREAM DISTANCE 


pi FOR lO-VORTEX CALCULATION OF FIGURE 3 
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.47781 

.45491 

.42444 

.38770 

.34598 

.30053 

.25259 

4.70670 

4.70670 

4.72001 

4.75109 

4.79233 

4.83678 

4.87955 

4.91594 

4.94209 

4.95581 

4.95656 

4.94484 

4.92163 

4.88814 

4.84576 

4.79606 

4.74074 

4.68156 

4.62022 

4.55827 

4.49707 

0 

-.03563 

-.06016 

-.06926 

-.06294 

-.04277 

-.00890 

.03732 

.09296 

.15421 

.21742 

.27962 

.33851 

.39225 

.43934 

.47858 

.50910 

.53041 

.54240 

.54528 

.53950 

4.50530 

4.50530 

4.51151 

4.52619 

4.54893 

4.57839 

4.61361 

4.65385 

4.69812 

4.74471 

4.79101 

4.83389 

4.87047 

4.89880 

4.91792 

4.92760 

4.92793 

4.91918 

4.90174 

4.87614 

4.84306 

0 

-.03633 

-.06944 

-. 0976 c 

-.11897 

-.13322 

-.14027 

-.13972 

-.13073 

-.11227 

-.08360 

-.04492 

.00237 

.05615 

.11415 

.17440 

.23525 

.29531 

.35331 

.40808 

.45858 

4.24130 

4.24130 

4.24436 

4.25135 

4.26259 

4.27818 

4.29810 

4.32229 

4.35062 

4.38291 

4.41894 

4.45842 

4.50097 

4.54606 

4.59298 

4.64079 

4.68832 

4.73428 

4.77735 

4.81635 

4.85032 

0 

-.03728 

-.07355 

-.10814 

-.14019 

-.16931 

-.19549 

-.21884 

-.23940 

-.25704 

-.27141 

-.28199 

-.28813 

-.28912 

-.28426 

- .27297 
-.25487 
-.22990 
-.19832 
-.16071 
-.11785 

3.90130 

3.90130 

3.90287 

3.90634 

3.91191 

3.91972 

3.92984 

3.94224 

3.95683 

3.97350 

3.99217 

4.01282 

4.03544 

4.06005 

4.08668 

4.U534 

4.14601 

4.17867 

4.21324 

4.24963 

4.28769 

0 

-.03797 

-.07558 

-.11259 

-.14866 

-.18359 

-.21733 

-.24986 

-.28113 

-.31103 

-.33946 

-.36632 

-.39156 

-.41516 

-.43713 

-.45749 

-.47623 

-.49333 

-.50870 

-.52223 

-.53378 

3.45730 

3.45730 

3.45808 

3.45976 

3.46244 

3.46618 

3.47102 

3.47696 

3.48398 

3.49207 

3.50123 

‘ 3.51146 

3.52277 

3.53514 

3.54855 

3.56297 

3.57835 

3.59466 

3.61186 

3.62993 

3.64884 

0 

-.03784 

-.07555 

-.11303 

-.15015 

-.18683 

-.22301 

-.25866 

-.29372 

-.32811 

-.36179 

-.39472 

-.42688 

-.45828 

-.48893 

-.51883 

-.54797 

-.57633 

-.60389 

-.63062 

-.65649 

2.56400 

2.86400 

2.86441 

2.86527 

2.86662 

2.86848 

2.87088 

2.87381 

2.87728 

2.88129 

2.88586 

2.89099 

2.89669 

2.90296 

2.90980 

2.91719 

2.92513 

2.93361 

2.94262 

2.95216 

2.96222 

0 

-.03765 

-.07525 

-.11276 

-.15012 

-.18732 

-.22432 

-.26UO 

-.29762 

-.33384 

-.36975 

-.40533 

-.44055 

-.47541 

-•50991 

-.54402 

-.57774 

-.61105 

-.64392 

-.67637 

-.70838 

1.95200 

1.95200 

1.95219 

1.95258 

1.95319 

1.95402 

1.95509 

1.95639 

1.95792 

1.95969 

1.96171 

1.96398 

1.96650 

1.96928 

1.97232 

1.97562 

1.97918 

1.98300 

1.98708 

1.99143 

1.99604 

0 

-.03526 

-.07050 

-.10571 

-.14086 

-.17595 

-.21096 

-.24587 

-.28067 

-.31535 

-.34990 

-.38431 

-.41856 

-.45265 

-.48658 

-.52033 

-.55390 

-.58728 

-.62046 

-.65343 

-.68619 


TABLE II.- VORTEX COORDINATES A, VERSUS DOWNSTREAM DISTANCE pi FOR lO-VORTEX CALCULATION OF FIGURE 4 



\ 



^2 

^3 

^3 



^5 






^8 

^8 


TI9 

Ajjo 

’lio 

0 

.5 

1.0 

1.5 
2.0 

2.5 

3.0 

3.5 

4.0 

4.5 

5.0 

5.5 

6.0 

6.5 

7.0 

7.5 

8.0 

8.5 
9.0 

9.5 
10.0 

4.89330 

4.89330 

4.89174 

4.88861 

4.88392 

4.87766 

4.86981 

4.86035 

4.84923 

4.83642 

4.82185 

4.80547 

4.78719 

4.76692 

4.74457 

4.72004 

4.69325 

4.66417 

4.63289 

4.59970 

4.56522 

0 

.03242 

.06476 

.09698 

.12903 

.16088 

.19246 

.22373 

.25462 

.28505 

.31495 

.34419 

.37265 

.40018 

.42655 

.45149 

.47464 

.49553 

.51358 

.52808 

.53831 

4.67200 

4.67200 

4.67160 

4.67080 

4.66961 

U . 66804 

4.66612 

4.66387 

4.66135 

4.65860 

4.65569 

4.65270 

4.64971 

4.64684 

4.64420 

4.64192 

4.64011 

4.63883 

4.63803 

4.63744 

4.63645 

0 

.01934 

.03868 

.05800 

.07729 

.09652 

.11568 

.13475 

.15374 

.17264 

.19147 

.21025 

.22902 

.24786 

.26687 

.28622 

.30614 

.32696 

.34914 

.37318 

.39962 

4.45870 

4.45870 

4.45^7 

4.45862 

4.45854 

4.45844 

4.45833 

4.45822 

4.45813 
4.45807 
4.45807 

4.45814 
4.45830 
4.45857 
4.45896 

4.45948 

4.46013 

4.46091 

4.46180 

4.46280 

4.46391 

0 

.01073 

.02147 

.03222 

.04296 

.05369 

.06441 

.07511 

.08580 

.09^7 

.10714 

.11781 

.12850 

.13924 

.15005 

.16097 

.17205 

.18334 

.19491 

.20685 

.21924 

4.24400 

4.24400 

4.24409 

4.24428 

4.24456 

4.24494 

4.24543 

4.24605 

4.24681 

4.24772 

4.24881 

4.25011 

4.25164 

4.25344 

4.25555 

4.25800 

4.26083 

4.26409 

4.26783 

4.27210 

4.27694 

0 

.00422 

.00845 

.01269 

.01693 

.02117 

.02540 

.02962 

.03383 

.03802 

.04220 

.04635 

.05049 

.0^61 

.05873 

.06284 

.06695 

.07109 

.07526 

.07951 

.08386 

02530 

4.02530 

4.02530 

4.02529 

4.02526 

4.02521 

4.02512 

4.02497 

4.02474 

4.02440 

4.02391 

4.02324 

4.02234 

4.02117 

4.01968 

4.01782 

4.01553 

4.01276 

4.00946 

4.00559 

4.00110 

0 

-.00095 

-.00190 

-.00283 

-.00373 

-.00458 

-.00536 

-.00606 

-.00666 

-.00715 

-. 00^53 

-.00778 

-.00790 

-.00788 

-.00772 

-.00742 

-.00699 

-.00644 

-.00578 

-.00505 

-.00427 

3.80930 

3.80930 

3.80944 

3.80971 

3.81011 

3.81063 

3.81125 

3.81196 

3.81273 

3.81353 

3.81433 

3.81510 

3.81580 

3.81639 

3.81684 

3.81711 

3.81717 

3.81700 

3.81658 

3.81590 

3.81497 

0 

-.00712 

-.01424 

-.02134 

-.02840 

-.03541 

-.04234 

-.04918 

-.05592 

-.06256 

-.06908 

-.07549 

-.08179 

-.08799 

-.09410 

-.10015 

-.10617 

-.11220 

-.11829 

-.12451 

-.13094 

3.59200 

3.59200 

3.59250 

3.59349 

3.59498 

3.59694 

3.59936 

3.60222 

3.60549 

3.60912 

3.61308 

3.61731 

3.62174 

3.62630 

3.63091 

3.63546 

3.63984 

3.64393 

3.64761 

3.65079 

3.65345 

0 

-.01471 

-.02941 

-.04408 

-.05869 

-.07322 

-.08766 

-.10198 

-.11617 

-.13021 

-.14411 

-.15787 

-.17149 

-.18498 

-.19839 

-.21176 

-.22519 

-.23882 

-.25283 

-.26751 

-.28320 

3.36000 

3.36000 

3.36128 

3.36385 

3.36771 

3.37287 

3.37935 

3.38718 

3.39639 

3.40703 

3.41914 

3.43279 

3.44805 

3.46500 

3.48373 

3.50434 

3.52695 

3.55166 

3.57857 

3.60775 

3.63918 

0 

-.02469 

-.04933 

-.07393 

-.09848 

-.12297 

-.14741 
-.17176 
-.19602 
- .22015 

-.24412 

-.26787 

-.29133 

-.31443 

-.33704 

-.35900 

-.38013 

-.40017 

-.41878 

-.43553 

-.44988 

2.85600 

2.85600 

2.85618 

2.85654 

2.85709 

2.85783 

2.85875 

2.85987 

2.86118 

2.86269 

2.86440 

2.86631 

2.86843 

2.87076 

2.87331 

2.87608 

2.87908 

2.88231 

2.88578 

2.88949 

2.89344 

0 

-.01778 

-.03555 

-.05329 

-.07099 

-.08863 

-.10618 

-.12363 

-.14097 

-.15817 

-.17522 

-.19212 

-.20884 

-.22538 

-.24174 

-.25791 

-.27388 

-.28965 

-.30523 

-.32062 

-.33583 

1.94400 

1.94400 

1 . 944 i 4 

1.94442 

1.94484 

1.94539 

1.94608 

1.94690 

1.94785 

1.94892 

1.95011 

1.95141 

1.95282 

1.95433 

1.95594 

1.95764 

1.95942 

1.96128 

1.96322 

1.96523 

1.96730 

0 

-.00967 

-.01933 

-.02899 

-.03861 

-.04819 

-.05771 

-.06717 

-.07656 

-.08586 

-.09507 

-.10417 

-.11317 

-.12205 

-.13081 

-.13944 

-.14793 

-.15628 

-.16450 

-.17257 

-.18050 
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TABLE III.- SINGLE WING VORTEX FORM 
(a) a = 5°; 9 = 0°; M = 2; r = 3/L; ii = 0; An = 1 


■p- 

co 


Vortex-induced velocities 


Vortex 
No. j 

rj 

Vo 


’'j 

H H 
1 1 

©.r 

’ll - ’ij 

’ll -® 


^1 

@x®/@ 

h 

-@x®/© 

2«r 



© 

@ 

® 


© 

® 

® 

@ 


© 

1 

0.12796 

1.50870 

0 

0 

0 

0 

0 

0 



2 

-.12796 

.66284 

0 

.84586 

0 

.71548 

-.15128 

0 



3 

.12796 

-.66284 

0 

2.17154 

0 

4.71559 

.05893 

0 

-0.02860 

0 

k 

- . 12796 

-1.50870 

0 

3.01740 

0 

9.10470 

-.04241 

0 




Crossflov/ velocities 


+ fi tan a 

^1^ 


(® ^ 

© - ® 

- 2 \ X @ 

a X @/@ 

a X ©/© 

I> 

An(@+@) 

’ll + 0 


© 


© 

© 

© 

© 

© 

© 

0 

0 

0 

0 

0 

2.27618 

0 

5.18100 

2.27618 

0 

0.03834 

0 

0.00974 

0 

0.00974 

1.50870 


Image vortex position 


U + A^I 

tan a 

@ X 0 

x^+ 0 

’ll + ® 

0 + 0 


^2 = ©/© 

©/© 

’I2 = - ©-^® 


0 

0 

0 


0 


© 

0 


1 

0.08749 

0.08749 

1.50870 

0.00974 

0.09723 

2.28563 

0.66008 

0.04254 

-0.04495 
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TABLE III.- SINGLE WING VORTEX FORM 
("b) a = 5°; T = 0°; M = 2; r = 3 A; = 1; An =1- Concluded 



Vortex induced velocities 

Vortex 
No. J 

Vo 


’'J 

\ 

>^1 - (3) 

’ll - ’■j 
-® 


'll 

©X ©/© 

-©X©/® 



© 

© 

@ 

© 

© 

© 


© 




1 

0.12796 

1.50870 

0.00974 

0 

0 

0 

0 

0 



2 

-.12796 

.66008 

-.04495 

.84862 

.05469 

•72315 

-.15016 

.00968 



3 

.12796 

-.66008 

-.04495 

2.16878 

.05469 

4.70660 

.05896 

-.00149 

-0.02835 

-0.00174 

4 

-.12796 

-1.50870 

.00974 

3.01740 

0 

9-10470 

-.042lu 

0 




Crossflow velocities 


+ P tan a 


& 

(© + @7 


-2\ X @ 

a X @/@ 

a X @/@ 



'll + © 

+ @ 



© 

© 

© 

© 




© 

0 

© 

0.09723 

2.27618 

0.00945 

5.22410 

2.26673 

-0.029338 

0.03787 

-0.00500 

0.00952 

-0.00674 

0.01926 

1.50196 


Image vortex position 


Ap 


tan a 


0.087^9 


@ 


>>2 = © / ® 


= 


0.17498 


1.50196 


0.01926 


0.19424 


2.29361 


0.65485 


0.08469 


-0.09029 
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(a) Vapor-screen photograph. 


A-18224 



Figure 1.- Vortex wake behind wing of high-aspect-ratio supersonic wing 
and body combination at angle of attack. 
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Figure 2.- Geometry of the flov field. 
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Aspect ratio 


2/3 



Figure 3 .- Rolling up of a vortex sheet "behind wing of subsonic leading-edge wing and body- 

configuration. 
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Figure 4 .- Rolling up of a vortex sheet behind wing of supersonic leading-edge wing and body 

conf iguration. 
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Figure 5 *- Rolling up of a vortex sheet segment containing eight equal-strength, equally spaced 

vortices . 
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0 12 3 4 

Span wise coordinate, A 



(a) Aspect ratio = 2/3 

Figure 6 .- Contribution of wing vortices, image vortices, and body cross- 
flow to total downvrash at tail location for various aspect ratios. 
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("b) Aspect ratio = 2 
Figure 6.- Continued. 
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Spanwise coordinate, A 



(c) Aspect ratio = 4 
Figure 6 .- Concluded. 


Spanwise coordinate, Spanwise coordinate 
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Figure 


"X 



//L ^ 1 1 1 1 L 

0 2 4 6 8 to !2 


Streamwise coordinate, /i 

(a) Aspect ratio = 2/3, M = 2 , r/sm = 0.6 




. - Vortex paths hehind two wing -body combinations at various 
angles of attack. 


Vertical coordinate. 


Spanwise coordinate, A Vertical coordinate 
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(c) Replot of figures 7 (a) and (b) against pa. 


Figure 7 -- Concluded. 
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Spanwise coordinate, X 



(a) a = 5 ° 


Figure 8 .- Contribution of wing vortices, image vortices, and body cross- 
flow to total downwash at tail location. 
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Spanwise coordinate^ A 



(Id) a = 10° 


Figure 8 .- Continued. 
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Semispan coordinate, A 


(c) a = 15° 


Figure 8 .- Continued. 
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Figiire 8 .- Continued. 
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Spanwise coordinate^ A 



(e) a = 11 ° 


Figure 8 .- Continued. 
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Span wise coordinate , A 
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(f) a = 16-1/2° 


Figure 8 .- Concluded. 
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Figure 9*“ Vortex paths behind wings of 2/3-aspect-ratio triangular -wing 
and. body conibinations with various ratios of body radius to wing 
semi span r/sm. 
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Figure 10.- Variation of single -vortex strength with body radius to wing 
semispan ratio for subsonic leading-edge configuration. 
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Figure 11 .- Contribution of wing vortices, image vortices, and body cross- 
flow to total downwash at tail location. 
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Figure 11.- Concluded. 
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-body 


Figure 12.- Comparison of rolling up of vortex sheet behind wing -body- 
combination (A = 2/3, r/sm = 0.2) with wing alone. 
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Figure 13* - Illustration of coordinates and dimensions of general wing- 

body configuration. 
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Figure l 4 .- Spanvise position of single vortex approximating circulation distribution. 
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Figure 15 - - Downvash along line y = z = 0 behind triangular wings, from reference ?• 
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Figure 20.- Approximated (spiral) and exact (circiaar) paths of tvo vortices of unequal strength. 
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Figure 21.- Vortex position error versus Ap at stations p = 1, 2, 3, 

and 4 for single wing vortex. 
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(a) Stations M- = 1 and 2. 

Figure 22.- Vortex position error versus A|i at stations |i = 1 , 2, 3 
and 4 for three wing vortex representations. 


Spanwise coordinate X, and lOxl^/^j 
(ij, vertical coordinate) 
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(■fa) Stations M’ = 3 and 4. 


Figure 22.- Continued. 
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Spanwise coordinate, A 


(c) Vortex paths from p = 0 to p = 4, with Ap = l/k, l/2, 1, and 2, viewed in crossflow plane. 


Figure 22.- Concluded. 
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Figure 23.- Comparison of stepwise approximation with Villat (exact 
crossflow) solution for two wing-body combinations (ref. 18) . 
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(a) Aspect ratio = 2/3, r/sj^ = 0.6 

Figure 24.- Effect of number of vortices on downwash at tail location 

for missile-type configuration. 
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Spanwise coordinate, X 



(b) Aspect ratio = k, r/sjj^ = 0.2 
Figure 2k.- Concluded. 
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(a) Aspect ratio = 2/3, = 0.6 

Figure 25.- Effect of niomber of vortices on downwash at tail location 

for airplane-type configuration. 
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(Td) Aspect ratio = 4, r/sjjj = 0.2 
Figure 25.- Concluded. 
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Aspect-ratio correction factor k. 
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Figure 27 -- Forecone of integration for equation (b 8 ) . 
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Figure 28.- Comparison of corrected slender-tody and linear theories at Mach line on wing 

trailing edge. 
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